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We study instabilities occurring in the electron system whose Fermi surface has flat regions on 
its opposite sides. Such a Fermi surface resembles Fermi surfaces of some high-Tc superconduc- 
tors. In the framework of the parquet approximation, we classify possible instabilities and derive 
renormalization-group equations that determine the evolution of corresponding susceptibilities with 
decreasing temperature. Numerical solutions of the parquet equations are found to be in qualitative 
agreement with a ladder approximation. For the repulsive Hubbard interaction, the antiferromag- 
netic (spin-density-wave) instability dominates, but when the Fermi surface is not perfectly flat, the 
d-wave superconducting instability takes over. 
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I. INTRODUCTION 



An interacting electron gas in one dimension has many unusual properties, such as the spin-charge separation, the 
I ' power law of correlation functions, and the linear dependence of the electron relaxation rate on temperature and 
, frequency (see Ref. for a review). These one-dimensional (ID) results are well established, in many cases exactly, 
' by applying a variety of mathematical methods including the Bethe Ansatz, the bosonization, and the parquet, or 
\^ ' the renormalization group. To distinguish the exotic behavior of the ID electron gas from a conventional Fermi-liquid 
On behavior, Haldane introduced a concept of the so-called Luttinger liquid ||]. 

The discovery of high-Tc superconductivity renewed interest in the Luttinger-liquid concept. Anderson suggested 
that a two-dimensional (2D) electron gas behaves like the ID Luttinger liquid, rather than a conventional Fermi 
S liquid . It is difficult to verify this claim rigorously, because the methods that prove the existence of the Luttinger 
liquid in ID cannot be applied directly to higher dimensions. The Bethe Ansatz construction does not work in higher 
(— ( dimensions. The bosonization in higher dimensions P-pT| converts a system of interacting electrons into a set of 
Q ' harmonic oscillators representing the electron density modes. This procedure replaces the exact Woo commutation 
O relations |^] with approximate boson commutators, which is a questionable, uncontrolled approximation. On the 
^ other hand, the parquet method, although not being as exact as the two other methods, has the advantage of being 
• ^ formulated as a certain selection rule within a standard many-body diagram technique; thus, it can be applied to 
higher dimensions rather straightforwardly. The parquet method has much in common with the renormalization-group 
^ , treatment of Fermi liquids [|2| . 
■ - - ' The ID electron gas has two types of potential instabilities: the superconducting and the density-wave, which 
manifest themselves through logarithmic divergences of the corresponding one-loop susceptibilities with decreasing 
temperature. Within the parquet approach, a sum of an infinite series of diagrams, obtained by adding and inserting 
the two basic one-loop diagrams into each other, is calculated by solving a system of nonlinear differential equations, 
which are nothing but the renormalization-group equations p^ . This procedure was developed for the first time 
for meson scattering ||l^ and later was successfully applied to the ID electron gas [ p"5|Jl^ , as well as to the Kondo 
problem and the X-ray absorption edge problem ||l8|] . By considering both the superconducting and the density- 
wave instabilities on equal footing and adequately treating their competition, the parquet approximation differs from a 
conventional ladder (or mean-field) approximation, commonly applied in higher dimensions, where only one instability 
is taken into account. Under certain conditions in the ID case, the superconducting and density-wave instabilities 
may cancel each other, giving rise to a non-trivial metallic ground state at zero temperature, namely the Luttinger 
liquid. In this case, the parquet derivation shows that the electron correlation functions have a power-law structure, 
which is one of the characteristic properties of the Luttinger liquid [ p^|Jl9|] . One may conclude that the competition 
between the superconducting and density-wave instabilities is an important ingredient of the Luttinger liquid theory. 

In a generic higher-dimensional case, where density-wave instability does not exist or does not couple to supercon- 
ducting instability because of corrugation of the Fermi surface, the parquet approach is not relevant. Nevertheless, 
there are a number of higher-dimensional models where the parquet is applicable and produces nontrivial results. 
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These include the models of multiple chains without single-electron hopping ||2^] and with single-electron hopping 
but in a magnetic field as well as the model of an isotropic electron gas in a strong magnetic field |^,^. In 
all of these models, the electron dispersion law is ID, which permits to apply the parquet method; at the same time, 
the interaction between electrons is higher-dimensional, which makes a nontrivial difference from the purely ID case. 
The particular version of the parquet method used in these cases is sometimes called the "fast" parquet, because, in 
addition to a "slow" , renormalization-group variable, the parquet equations acquire supplementary, "fast" variables, 
which label multiple electron states of the same energy. 

Taking into account these considerations, it seems natural to start exploring a possibility of the Luttinger liquid 
behavior in higher dimensions by considering a model that combines ID and higher-dimensional features. This is 
the model of an electron gas whose Fermi surface has flat regions on its opposite sides. The flatness means that 
within these regions the electron dispersion law is ID: The electron energy depends only on the one component of 
momentum that is normal to the flat section. On the other hand, the size of the flat regions is finite, and that property 
differentiates the model from a purely fD model, where the size is infinite, since nothing depends on the momenta 
perpendicular to the direction of a ID chain. A particular case of the considered model is one where the 2D Fermi 
surface has a square shape. This model describes 2D electrons on a square lattice with the nearest-neighbor hopping 
at the half filling. It is a simplest model of the high-Tc superconductors. 

The model has already attracted the attention of theorists. Virosztek and Ruvalds studied the "nested Fermi liquid" 
problem within a ladder or mean- field approximation [ p4U25| ]. Taking into account the ID experience, this approach 
may be considered questionable, because it does not treat properly the competition between the superconducting 
and the density-wave channels. Houghton and Marston mapped the flat parts of the Fermi surface onto discrete 
points. Such an oversimplification makes all scattering processes within the flat portion equivalent and artificially 
enhances the electron interaction. Mattis and Hlubina used the bosonization to treat the interaction between 
the electron density modes and claimed to solve the model exactly. However, mapping of the flat Fermi surface onto 



I indicated that the treatment of Mattis and Hlubina 
app scattering on different quantum chains require a 



quantum chains and subsequent bosonization by Luther [£8 
is insufficient, because the operators of backward and umk 
consistent renormalization-group treatment. Luther did not give solution to this problems, as well as he missed the 
interaction between the electrons located on four different quantum chains. 

In the present paper, we solve the model consistently, using the fast parquet approach, where all possible instabilities 
occurring in the electron system with the flat regions on the Fermi surface are treated simultaneously. This approach 
was applied to the problem earlier ]29|] in order to explain the antiferromagnetism of chromium. In the present paper, 
we advance the study further by including the order parameters of the odd symmetry, missed in performing 
detailed numerical calculations, and investigating the effect of a curvature of the Fermi surface. To simplify numerical 
calculations and to relate to the high- Tc superconductors, we consider the 2D case, although the method can be 
straightforwardly generalized to higher dimensions as well. 

We find that the presence of the boundaries of the flat portions of the Fermi surface has a dramatic effect on 
the solutions of the parquet equations. Even if the initial vertex of interaction between electrons does not depend 
on the momenta along the Fermi surface (which are the "fast" variables), the vertex acquires a strong dependence 
on these variables upon renormalization, which greatly reduces the feedback coupling between the superconducting 
and density-wave channels relative to the ID case. Instead of the two channels canceling each other, the leading 
channel, which is the spin-density-wave (SDW) in the case of the repulsive Hubbard interaction, develops its own 
phase transition, inducing on the way a considerable growth of the superconducting d-wave susceptibility. At the same 
time, the feedback from the superconducting to the SDW channel, very essential in the ID case, is found negligible 
in the 2D case. These results are in qualitative agreement with the picture of the antiferromagnetically-induced 
d-wave superconductivity, which was developed within a ladder approximation for the flat Fermi surface in p5| and 
for a generic nested Hubbard model in | |30[ |. Recent experiments strongly suggest that the high- Tc superconductivity 
is indeed of the d-wave type On the other hand, our results disagree with Refs. [ p6|j27| ]. The origin of the 

discrepancy is that the bosonization arbitrarily replaces the exact Woo commutation relations |^ by approximate 
boson commutators; thus, the renormalization of the electron-electron interaction, which is an important part of the 
problem, becomes neglected. 

In addition to having the flat sides, the square Fermi surface also has sharp corners, where the saddle points of the 
electron dispersion law, which produce the van Hove singularity in the density of states, are located. The presence of 
the van Hove singularity at the Fermi level enhances the divergence of the superconducting and density-wave loops 
to the square of the temperature logarithm. The fast parquet problem was formulated in this case in Ref. |Q, where 
the contribution from the flat sides, being less divergent than the contribution from the saddle points, was neglected. 
The present paper completes the study by considering a Fermi surface with the flat sides and rounded corners, that 
is, without saddle points at the Fermi level. Our physical conclusions for both models are in qualitative agreement. 

As photoemission experiments demonstrate (see also ) , many of the high- Tc superconductors indeed have 
flat regions on their Fermi surfaces. Hence, some of the results of this paper may be applicable to these materials. 
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However, the primary goal of our study is to elucidate general theoretical concepts rather than to achieve detailed 
description of real materials. 

In order to distinguish the new features brought into the problems by introducing higher dimensions, we present 
material in an inductive manner. In Sec. we recall the derivation of the parquet equations in the simplest case of 
ID spinless electrons. In Sec. [II, we generalize the procedure to the case of ID electrons with spin |T5|JT^ ]. Then, we 
derive the parquet equations in the 2D case in Sec. IV and solve them numerically in Sec. 0. The paper ends with 
conclusions in Sec. VI. 



II. PARQUET EQUATIONS FOR ONE-DIMENSIONAL SPINLESS FERMIONS 

Let us consider a ID electron gas with a Fermi energy /x and a generic dispersion law e{kx), where e is the energy 
and kx is the momentum of the electrons. As shown in Fig. |l|, the Fermi surface of this system consists of two points 
located at kx = ikp, where hp is the Fermi momentum. Assuming that the two points are well separated, let us 
treat the electrons whose momenta are close to ikp as two independent species and label them with the index ±. In 
the vicinity of the Fermi energy, the dispersion laws of these electrons can be linearized: 

e±{kx) ^ ±vpkx, (2.1) 

where the momenta kx are counted from the respective Fermi points ikp for the two species of the electrons, ±vf 
are the corresponding Fermi velocities, and the energy e is counted from the chemical potential /i. 

First, let us consider the simplest case of electrons without spin. The bare Hamiltonian of the interaction between 
the ± electrons, i/inti can be written as 

ij,, = 9 1 ^^^^tiki'^ + ki'^ - k^')i^tmi^^mi^+m. (2.2) 

where g is the bare vertex of interaction, and the operators V-'i and '4)± create and destroy the ± electrons. 

The tendencies toward the superconducting or density-wave (2^^) instabilities in the system are reflected by the 
logarithmic divergences of the two one-loop diagrams shown in Fig. y, where the solid and dashed lines represent the 
Green functions G+ and of the + and — electrons, respectively. The two diagrams in Fig. |^ differ in the mutual 
orientation of the arrows in the loops. In the Matsubara technique, the integration of the Green functions over the 
internal momentum kx and energy u;„ produces the following expressions for the two diagrams: 



/dk 
'^x ~^ Qx) 

n 



dk^r 1 



27r {iuJn + VFkx){iuJn + i^m - VF{kx + qx)) 



^ Inf „ , ,^ n (2.3) 



2t:vf \ max{T, \vFqx \,\^m\} 

where the upper sign corresponds to the superconducting and the lower to the density-wave susceptibility. In Eq. 
( ^.3D , T is the temperature, Vim is the external energy passing through the loop, and qx is the external momentum 
for the superconducting loop and the deviation from ikp for the density- wave loop. With logarithmic accuracy, the 



value of the integral (|2.3|) is determined by the upper and lower cutoffs of the logarithmic divergence. In Eq. (2.3), the 



upper and lower cutoffs are written approximately, up to numerical c oeffi cients of the order of unity, whose logarithms 



are small compared to ^ 3> 1. The variable ^, introduced by Eq. (2^), plays a very important role in the paper. 
Since ^ is the logarithm of the infrared cutoff, the increase of ^ represents renormalization toward low temperature 
and energy. 

The two primary diagrams of Fig. ^ generate higher-order corrections to the vertex of interaction between electrons, 
7, as illustrated in Fig. |^. In this Figure, the dots represent the bare interaction vertex g, whereas the renormalized 
vertex tIs shown as a circle. The one-loop diagrams in Fig. ^ are the same as in Fig. ^j. The first two two-loop diagrams 
in Fig. |3| are obtained by repeating the same loop twice in a ladder manner. The last two diagrams are obtained 
by inserting one loop into the other and represent coupling between the two channels. The diagrams obtained by 
repeatedly adding and inserting the two basic diagrams of Fig. |^ in all possible ways are called the parquet diagrams. 
The ladder diagrams, where only the addition, but not the insertion of the loops is allowed, represent a subset of the 
more general set of the parquet diagrams. Selection of the parquet diagrams is justified, because, as one can check 
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calculating the diagrams in Fig. ^, they form a series with the expansion parameter g^: 7 = gJ2'^=o ^nigO""- If the 
bare interaction vertex g is small and the temperature is sufhciently low, so that ^(T) is big, one can argue jlj-jl^] 
that nonparquet diagrams may be neglected, because their expansion parameter g is small compared to the parquet 
expansion parameter g^. 

Every diagram in Fig. |3| except the bare vertex g, can be divided into two disconnected pieces by cutting one solid 
and one dashed line, the arrows of the cut lines being either parallel or antiparallel. The sum of those diagrams where 
the arrows of the cut lines are parallel (antiparallel) is called the superconducting (density-wave) "brick". Thus, the 
vertex 7 can be decomposed into the bare vertex g, the superconducting brick C, and the density-wave brick Z: 



j^g + C + Z. 



(2.4) 



Eq. (^^) is illustrated in Fig. ^, where the bricks are represented as rectangles whose long sides, one being a solid 
and another a dashed line, represent the lines to be cut. 



In a general case, the vertices and the bricks depend on the energies and momenta (uji, lu2, lu^, vpk 



(1) 



and 



vpkx ) of all incoming and outgoing electrons. Equations for the bricks can be found in closed form in the case where 
all their arguments are approximately equal within the logarithmic accuracy, that is, the ratios of the arguments and 
of their linear combinations are of the order of unity p^-p^ . Practically, this means that all vertices and bricks are 
considered to be functions of the single renormalization-group variable ^, defined in Eq. (2.3). It was proved in p4 ] 
that the two pieces obtained by cutting a brick are the full vertices of interaction, as illustrated graphically in FigTlf . 
Analytically, the equations for the bricks are 



rfC7(C)7(C), 



7(07(0- 



(2.5a) 
(2.5b) 



The two vertices 7 in the r.h.s. of Eqs. (2.5) represent the two pieces obtained from a brick by cutting, whereas 
the integrals over ^ represent the two connecting Green functions being integrated over the internal momentum and 
energy of the loop. The value of the renormalized vertex 7(0 changes as the integration over ^ progresses in Eqs. 
( |2.5[ ). In agreement with the standard rules of the diagram te chniq ue a pair of the parallel (antiparallel) lines in 
Fig. H produces a negative (positive) sign in the r.h.s. of Eq. ( 2.5a ) [( |2.5b| )]. 
Eqs. (|2.5[) can be rewritten in differential, renormalization-group form: 



dm 



= -7(07(0, 
= 7(07(0, 



C(e=0) = 0; 
Z(^=0) = 0. 



Combining Eqs. ( p.q ) with Eq. (2^), we find the renormalization equation for the full vertex 7: 

^7(0 



d^ 

7(^-0) 



7(07(0 - 7(07(0 = 0, 
= 9- 



(2.6a) 
(2.6b) 

(2.7a) 
(2.7b) 



We see that the two terms in the r.h.s. of Eq. ( ^.7a ), representing the tendencies toward density-wave and supercon- 
ducting instabilities, exactly cancel each other. In a ladder approximation, where only one term is kept in the r.h.s., 
the result would be quite different, because 7(0 would diverge at a finite ^ indicating an instability or generation of 
a pseudogap in the system. 

In order to study possible instabilities in the system, we need to calculate corresponding generalized susceptibilities. 
For that purpose, let us add to the Hamiltonian of the system two fictitious infinitesimal external fields h^c and /idw 
that create the electron-electron and electron-hole pairs: 



^^Gxt — 



dqx dkx 
2n 27r 



kr. 



+ /iDw(gx) '01 [kx + 



I) ('- 



qx 
2 



H.c. 



(2.8) 



Now we need to introduce triangular vertices Tsc and Tdvv that represent the response of the system to the fields 
/isc and /idw- Following the same procedure as in the derivation of the parquet equations for the bricks |15 l^, |2^ , p9t , 
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we find the parquet equations for the triangular vertices in graphic form, as shown in Fig. ^. In that Figure, the 
filled triangles represent the vertices Tsc and Tdwi whereas the dots represent the fields hsc and /idw- The circles, 
as in the other Figures, represent the interaction vertex 7. Analytically, these equations can be written as differential 
equations with given initial conditions: 



7(0^DW(0, 



^c(O) = hsc; 

7dw(0) = /iDW- 



(2.9a) 
(2.9b) 



We will often refer to the triangular vertices T as the "order parameters" . Indeed, they are the superconducting and 
density- wave order parameters induced in the system by the external fields /igc and /idw- for a finite hi (i=SC, 
DW), a vertex 7i(^), which is proportional to hi, diverges when ^ ^ ^c, this indicates that a spontaneous order 
parameter appears in the system, that is, the order parameter may have a finite value even when the external field 
hi is zero. The external fields are introduced here on ly as auxiliary tools and are equal to zero in real systems. We 
also note that the two terms in the r.h.s. of Eq. (2.8) are not Hermitially self-conjugate; thus, the fields hi are the 



complex fields. Consequentl y, th e order parameters are also complex, so, generally speaking, T and T* do not 
coincide. According to Eqs. (2/7), 7(^) = g, so Eqs. (2_^) have the following solution: 



TsciO = ^scexp(-50, 
7bw(0 = ^Dwexp(5r^). 



(2.10a) 
(2.10b) 



Now we can calculate the susceptibilities. The lowest order corrections to the free energy of the system due to the 
introduction of the fields hsc and /idWj ^SC and -Fbwj obey the parquet equations shown graphically in Fig. |7| and 
analytically below: 



FsciO 



dcrsc(C)rs*c(C), 



Jo 



(2.11a) 
(2.11b) 



Substituting expressions (2.10) into Eqs. (2.11) and dropping the squares of hsc and /idw, we find the susceptibilities: 



XSciO = - [ exp(-2gO - 1 ] /2.g, (2.12a) 
XDwiO = [ cxp(2ga - 1 ] /25- (2.12b) 

According to Eqs. ( ^.12| ), when the interaction between electrons is repulsive (attractive), that is, g is positive 
(negative), the density- wave (superconducting) susceptibility increases as temperature decreases (T — > and ^ — > 00): 

/ fi \ 
Xdw(sc)(0 oc exp(±2gO = tt^i , m 11 ■ (2-13) 



Susceptibilities ( 2.13| ) have power dependence on the temperature and energy, which is one of the characteristic 



properties of the Luttinger liquid. The susceptibilities are finite at finite temperatures and diverge only at zero 
temperature, in agreement with the general theorem |]35|| that phase transitions are impossible at finite temperatures 
in ID system s. M athematically, the absence of divergence at finite ^ is due to the cancellation of the two terms in the 
r.h.s. of Eq. ( ^.7a| ) and subsequent nonrenormalization of 7(^). This nontrivial ID result can be obtained only within 
the parquet, but not the ladder approximation. 



III. PARQUET EQUATIONS FOR ONE-DIMENSIONAL FERMIONS WITH SPIN 

Now let us consider ID electrons with spin. In this case, there are three vertices of interaction, conventionally 
denoted as 71, 72, and 73, which represent backward, forward, and umklapp scattering, respectively [ [l5| , |l6[ . Umklapp 
scattering should be considered only when the change of the total momentum of the electrons in the interaction 
process, 4fcF, is equal to the crystal lattice wave vector, which may or may not be the case in a particular model. In 
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this paper, we do not consider the vertex 74, which describes the interaction between the electrons of the same type 
(+ or — ), because this vertex does not have logarithmic corrections. The bare Hamiltonian of the interaction, Hint, 
can be written as 



Hint = 



dJxix dkx dkx 



2tt 2tt 2tt 



(3.1) 



where the coefficients 51-3 denote the bare (unrenormalized) values of the interaction vertices 71-3. The operators 
and il^as create and destroy electrons of the type s = ± and the spin c = ti- The spin structure of the interaction 
Hamiltonian is dictated by conservation of spin. We picture the interaction vertices in Fig. |^, where the solid and 
dashed lines represent the + and — electrons. The thin solid lines inside the circles indicate how spin is conserved: 
The spins of the inc oming and outgoing electrons connected by a thin line are the same. According to the structure 
of Hamiltonian (3.1), the umklapp vertex 73 describes the process where two + electrons come in and two - electrons 
come out, whereas the complex conjugate vertex 73 describes the reversed process. 
The three vertices of interaction contain six bricks, as shown schematically in Fig. |^: 



71 = 51 + C*! + Zl, 

12 = 92 + C2 + Z2, 

13 = 93 + Zl + Zii, 



(3.2a) 
(3.2b) 
(3.2c) 



where Ci and C2 are the superconducting bricks, and Zi, Z2, Zj, and Zjj are the density- wave bricks. In Fig. ^, the 
thin solid lines inside the bricks represent spin conservation. The umklapp vertex has two density-wave bricks Zj and 
Zji, which differ in their spin structure. 

Parquet equations for the bricks are derived in the same manner as in Sec. ^ by adding appropriate spin structure 
dictated by spin conservation. It is convenient to derive the equations graphically by demanding that the thin spin 
lines are continuous, as shown in Fig. |l^. Corresponding analytic equations can be written using the following rules. 
A pair of parallel (antiparallel) lines connecting two vertices in Fig. |l^ produces the negative (positive) sign. A closed 
loop of the two connecting lines produces an additional factor —2 due to summation over the two spin orientations of 
the electrons. 

d^ 
dZiiO 

d^ 
dZ2{0 

dC 
dZiiO 

d( 
dZniO 

dC 



= -271(072(0, 


(3.3a) 


= -7?(0-7l(0, 


(3.3b) 


= 271(0 72(0-272(0, 


(3.3c) 


= 72(0+73(073(0, 


(3.3d) 


= 273(0172(0- 71 (0], 


(3.3e) 


= 273(0 72(0- 


(3.3f) 



Combining Eqs. (3.2) and (3.3), we obtain the well-known closed equations for renormalization of the vertices |16| |: 

d7i(0 



d^ 
dl2{0 

dS, 
dl3{0 

d^ 



= -27?(0, 
= -7?(0 + 73(073(0, 
= 273(0 [272 (0-71 (01- 



(3.4a) 
(3.4b) 
(3.4c) 
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In the presence of spin, the electron operators in Eq. ( p.8[ ) and, correspondingly, the fields hi and the triangular 
vertices 7i{^) acquire the spin indices. Thus, the superconducting triangular vertex TsciO becomes a vector: 



(3.5) 



Parquet equations for the triangular vertices are given by the diagrams shown in Fig. where the spin lines should 
be added in the same manner as in Fig. Wu. The superconducting vertex obeys the following equation: 



(3.6) 



where the matrix rsc(C) is 



rsc(e) 



/ -72 + 71 

-J/J 

V -72+71 



(3.7) 



Linear equation ( B.6 ) is diagonalized by introducing the singlet, Tssc j f^nd the triplet, Ttsc i superconducting triangular 
vertices: 



Tssc(C) = Tsc(C)-^sc(C), 



-Ih 



which obey the following equations: 

'^'^SC(TSC) (0 



[T7i(0 -72(0] ^SC(TSC)(0- 



(3.8a) 
(3.8b) 



(3.9) 



In Eq. ( |3.9|) the sign — and the index SSC correspond to the singlet superconductivity, whereas the sign + and the 
index TSC correspond to the triplet one. In the rest of the paper, we use the index SC where discussion applies to 
both SSC and TSC. 

Now let us consider the density- wave triangular vertices, first in the absence of umklapp. They form a vector 



(3.10) 



which obeys the equation 



with the matrix 



— ^DwiO %w(C) 



(3.11) 



DW 



( -71 + 72 -71 

72 

72 

\ -71 -71 + 72 



(3.12) 



Eq. ( |3.11 ) is diagonalized by introducing the charge-, TcdWi and the spin-, TsdW; density- wave triangular vertices: 
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TcuwiO = T^UO + ^Dw(0, (3.13a) 
TsDwiO - I r^liO 1 > (3-13b) 



which obey the following equations: 



[-271(0+ 72(0] ^CDw(0, (3.14a) 
72(0 %uw{0- (3.14b) 



When the iimklapp vertices 73 and 7I are introduced, they become ofFdiagonal matrix elements in Eqs. ( ^.14 ), 
mixing Tcdw and Tsdw with their complex conjugates. Assuming for simplicity that 73 is real, we find that the 
following linear combinations diagonalize the equations: 

^DW(SDW)± = ^DW(SDW) ± ^DW(SDW)' (3.15) 

and the equations become: 



d^ 

'^^DW±(0 

dC 



[-271(C) +72(0 T 73(C)] rcDW±(0, (3.16a) 
[72(0 ±73(0] ^SDW±(0- (3.16b) 



If the external fields hi are set to unity in the initial conditions of the type (2^) for all triangular vertices i ~ SSC, 
TSC, CDW +, an d SDW±, then the corresponding susceptibilities are equal numerically to the free energy corrections 
of the type ( ^.ll]) : 

X:{i)= t dCUOViC). (3.17) 



Eqs. (B.4), ( |3.9| ), ( 3.16| ), and (3.17) were solved analytically in Ref. [|T^], where a complete phase diagram of the ID 
electron gas with spin was obtained. 



IV. PARQUET EQUATIONS FOR TWO-DIMENSIONAL ELECTRONS 

Now let us consider a 2D electron gas with the Fermi surface shown schematically in Fig. ^ It contains two pairs 
of flat regions, shown as the thick lines and labeled by the letters a and h. Such a Fermi surface resembles the Fermi 



surfaces of some high-Tc superconductors |33|. In our consideration, we restrict the momenta of electrons to the flat 
sections only. In this way, we effectively neglect the rounded portions of the Fermi surface, which are not relevant 
for the parquet consideration, because the density-wave loop is not divergent there. One can check also that the 
contributions of the portions a and b do not mix with each other in the parquet manner, so they may be treated 
separately. For this reason, we will consider only the region a, where the 2D electron states are labeled by the two 
momenta kx and ky^ the latter momentum being restricted to th e in terval [—k'^\k'^^]. In our model, the energy 



of electrons depends only on the momentum k^ according to Eq. (2J) and does not depend on the momentum ky. 
We neglect possible dependence of the Fermi velocity Dp' on ky'^ it was argued in Ref. |2^ that this dependence is 
irrelevant in the renormalization-group sense. 

In the 2D case, each brick or vertex of interaction between electrons acquires extra variables k'if', and k^y^ in 

addition to the ID variables wi, uj2, W3, VFk^\ and VFk^\ These two sets of variables play very different 

roles. The Green functions, which connect the vertices and produce the logarithms depend only on the second set 
of variables. Thus, following the parquet approach outlined in the previous Sections, we dump all the uj and vpk^ 
variables of a vertex or a brick into a single variable At the same time, the ky^\ k^\ and ky^^ variables remain 
independent and play the role of indices labeling the vertices, somewhat similar to the spin indices. Thus, each vertex 
and brick is a function of several variables, which we will always write in the following order: "f{ky'^\ ^y^^i ^'y \ k'y^ ] 0- 
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It is implied that the first four variables satisfy the momentum conservation law k^y^ + k'^^ 



k^^ + k'i^' , and each of 

] . The assignment of the variables ky^' , k'y'' , k)^' , and k^^ to the ends of the 

vertices and bricks is shown in Fig. ^, where the labels kj {j — 1 — 4) should be considered now as the variables ky'^ . 
To shorten notation, it is convenient to combine these variable into a single four-component vector 

/C=(fcW,fc(^);fc(3),fcf), (4.1) 

so that the relation between the vertices and the bricks can be written as 



„(4) 



them belongs to the interval [— fcy°\ k'y^ 



(1) i,(2) i^d) 



7i(/C,e)=.9i+Ci(/C,0 + ^i(^,0, 
72(^,C)=.92 + C2(/C,e) + ^2(/C,0, 
73(/C,C)=53 + ^/(/C,e) + ^//(^,C)- 



(4.2a) 
(4.2b) 
(4.2c) 



After this introduction, we are in a position to write the parquet equations for the bricks. These equations are 
shown graphically in Fig. ^ where again the momenta kj should be understood as ky^ . Analytically, the equations 
are written below, with the terms in the same order as in Fig. pll: 



-7i(^i,C) °72(/C'i,C) - 72(^1,0 ° 71 (^'1,0 7 
-71 (^1 , ° 71 (/C'l , - 72 (/Ci , ° 72 (/C; , , 

71 (/C2 , ° 72 {IC', ,0+72 {IC2 , 6 ° 71 {IC'2 , - 271 (^2 , ° 71 , 

- 273 (/C2 , ° 73 (^2 : + 73 (/C2 , C) ° 73 (^2 , + 73 (/C2 , ° 73 (^2 , 6 , 

72 (^2 , ° 72 (/Cf, , a + 73 (/C2 , C) o 73 (/C^ , C) , 

73 (^3 , ° 72 (/C!, , C) + 72 (/C3 , ° 73 (/C^ , + 71 (^3 , ° 73 {IC', , 
+73 (/C3 , ° 71 (/C;,, - 273 (/C3 , C) o 71 itC'^ , ~ 271 (^3 , ° 73 ilC', , , 
73 (/C2 , ° 72 (/C^' , e) + 72 (/C2 , a ° 73 (/C^' , , 



where 



1^1 — yf^y J '^y T '^y I'^y J' '^1 " K'^y t ^y J '^y 7 '^y )i 

'^2 — \r<-y , l^y T '^y t'^v )' '^2 — y^y ' '^y ' '^y '""a )' '^2 — l"-?^ i "-j/ t '^y i "'y /: 

and the tilde and the bar operations are defined as 

^.('1,(1) 1.(2). 1,(3) 7,(4). f:\^^JL.(l} 1,(2). l(4) .(3). 

73(fc^"'"'': fcy^''; C) = 73 (^a'*^ ^^'^''i ^y^^^T^y^^'^ 0- 



(4.3a) 
(4.3b) 

(4.3c) 
(4.3d) 

(4.3e) 
(4.3f) 

(4.4a) 
(4.4b) 
(4.4c) 

(4.5a) 
(4.5b) 



In Eqs. (O), we introduced the operation o that represents the integration over the internal momenta of the loops 



in Fig. [1C|. It denotes the integration over the intermediate momentum ki/^^ with the restriction that both ky'^' and 
ky^\ another intermediate momentum determined by conservation of momentum, belong to the interval [—ky^\ k 
For example, the explicit form of the first term in the r.h.s. of Eq. (4.3a) is: 



(0) 1,(0)1 

■y J 



7i(^i,e)°72(/C'i,0 , 

<Ck^^^ -{-k^'^^ —k^^^ 'Ck^'^^ 

X7i(fc^^\ fc^^'; + \f - k^y^ 72(fc« + fcf - k^y^\ fc(-^); kf\kf- 0. 



(4.6) 



Eqs. ( |4.3| ) and (4.2) with definitions ( |4.l| ), (4.4), and (4.5) form a closed system of integrodifferential equations, which 
will be solved numerically in Sec. ^ The initial conditions for Eqs. (4.2) and ( |4.2| ) are that all the C and Z bricks 
are equal to zero at ^ = 0. 
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Parquet equations for th e superconducting triangular vertices can be found in the 2D case by adding the ky momenta 
to the ID equations (3_^). The equations are shown graphically in Fig. where the momenta k and q should be 



interpreted as ky and qy: 

d%SC(TSC){ky^%^^) r (V C\ n- 111 <:\ I A '7\ 
-Q-^ = JSSC(TSC)(/^SC,€) °^SSC(TSC)(S,%,?). (4-7) 

where 

/ssc(Tsc)(^sc,$) = T7i(^sc,0 -72(^sc,C): (4-8) 

^SC = (fc; + %/2, -k'y + qyl2- ~ky + qy/2, ky + 9^/2), (4.9) 

and the operator o denotes the integration over k'y with the restriction that both k'y + qy/2 and —ky+qy/2 belong to the 

interval [~kl^\ ^lul- ^^'^ T signs in front of 71 in Eq. ( |4.8| ) correspond to the singlet and triplet superconductivity. As 
discussed in Sec. ||, the triangular vertex 7sc(fcj,, %,£,) is the superconducting order parameter, qy and ky being the y- 
components of the total and the relative momenta of the electrons in a Cooper pair. Indeed, the vertex Tsc{ky, qy, S,) 
obeys the linear equation shown in Fig. ^, which is the linearized Gorkov equation for the superconducting order 
parameter. As the system approaches a phase transition, the vertex Tsc(ky, qy, S^) diverges in overall magnitude, 
but its dependence on ky for a fixed qy remains the same, up to a singular, ^-dependent factor. The dependence of 
Tsc{ky, qy,0 describes the distribution of the emerging order parameter over the Fermi surface. The numerically 

found behavior of T^ciky, qy, S,) is discussed in Sec. 0. 

Due to the particular shape of the Fermi surface, the vertices of interaction in our 2D model have two special 
symmetries: with respect to the sign change of all momenta ky and with respect to the exchange of the + and — 
electrons: 

7.(fc«,fc(^); fcf),fc(4); e) =7.(-4'^-fcf ; "^f '"^f 5 0, i = 1,2,3; (4.10a) 
7.(fc«,fc(^); 0=7.(fcf\fc^'^ ; 0, * = 1,2,3; (4.10b) 

73(fc«,fcf ; ; = 73(fcf 0, (4.iOc) 



where in Eq. (4.10c) we assume for simplicity that 73 is real. As a consequence of ( 4.10| ), Eqs. (4.7) are invariant with 



respect to the sign reversal of ky in TgQ{ky, qy, ^) at a fixed Qy. The following combinations of the triangular vertices 



form two irreducible representations of this symmetry, that is, they are independent and do not mix in Eqs. (4.7) 



'^SSC(TSC)('^a,9y,0 = ^SC(TSC)(^y,<?y,0 ± '?SSC(TSC) ("^y , <?y , 0- (4-11) 

The triangular vertices "Tsscctsc) (^2/, , ^) describe the superconducting order parameters that are either symmetric 
or antisymmetric with respect to the sign change of ky. When is extended over the whole 2D Fermi surface 

(see Fig. ^l]), it acquires the s-wave symmetry, whereas T^g^ the d-wave symmetry. The symmetrized vertices 



^^c(TSC)(^J" ^) ob^y the same Eqs. (4/7) as the unsymmetrized ones. 



The equations for the density-wave triangular vertices are obtained in a similar manner: 

/CDW±(^DW,6 °^CDW±(^y,9y,0, (4.12a) 
/SDW±(^DW,0 °^SDW±(^y'9a'C), (4.12b) 



where 



/CDW±(^DW,C) = -27l(^DW,C) T 273(/CdW,C) +72(^DW,0 ±73(^DW,C), (4-13) 
/SDW±(^DW,^) = 72(^DW,'?) ±73(^DW,C), (4-14) 
^DW = {k'y + qy/2, ky ~ qyl2; k'y - qy/2, ky + qy/2). (4.15) 



The ± signs in the subscripts of T in Eqs. (4.12) and in front of 73 in Eqs. (4.13)-(4.14) refer to the umklapp symmetry 
discussed in Sec. [II, whereas the ± signs in the superscripts of T refer to the symmetry with respect to sign reversal 
of ky, discussed above in the su per conducting case. The fc^-antisymmetric density waves are actually the waves of 
charge current and spin current [^,|2j , also known in the so-called flux phases ||3^ . 
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Once the triangular vertices 7^ are found, the corresponding susceptibiHties Xi calculated according to the 
following equation, similar to Eq. ( 3.17 ): 



dc 



dk. 



2tt 



(4.16) 



where the integrati on o v er ky is restr icted so that both ky ±qy/ 2 belong to the interva l \— ky , k 



(0) l(0)i 



Using functions ( |4.^ ), ( [4.13 ), and (4.14) and symmetries ( 1.1C| ), we can rewrite Eqs. (4.3) in a more compact form. 
For that purpose, we introduce the SSC, TSC, CDW, and SDW bricks that are the linear combinations of the original 
bricks: 



SSC(TSC) 



^CDW± — Z2 — 2Zi ± {Zji — 2Zi), 



where the tilde operation is defined in Eq. (4.5a). Then, Eqs. (4.3) become: 

5CSSC(TSC)(^,0 



dZsBw±{IC,0 



= ^/SSC(TSC)(^17 C) ° /sSC(TSC)(^'liOi 
= /CDW± (A^3 , ° /CDW± (A^3 , 7 
= /SDW±(^^2,C) °/sDW±(^2:C)- 



(4.17a) 

(4.17b) 
(4.17c) 



(4.18a) 
(4.18b) 
(4.18c) 



The parquet equations in the form ( 4.18 ) were obtained in Ref. ||2 



It is instructive to trace the difference between the parquet equations (4.18) a nd th e corresponding ladder equations. 
Suppose that, for some reason, only one brick, say Cssc: among the six bricks (|4.17| ) is appreciable, whereas the other 
bricks may be neglected. Using definitions (p^) and (4.8), we find that Eq. (4.18a) becomes a closed equation: 



g/ssc(/C,0 



= /ssc(^i,e)°/ssc(/C'i,0, 



where 



/ssc(A:i,C) = -.91 - 52 - C'ssc(/C,^). 



(4.19) 



(4.20) 



Eq. (4.19) is the lad der eq uation for the singlet superconductivity. When the initial value —{gi + g-i) of the vertex 
/ssc is positive, Eq. (4.19) has a singular solution (/ssc ^ 00 at ^ ^ ^c), which describes a phase transition into the 
singlet superconducting state at a finite temperature. Repeating this consideration for every channel, we construct 
the phase diagram of the system in the ladder approximation as a list of necessary conditions for the corresponding 
instabilities: 



SSC 
TSC 
CDW+ 
CDW- 
SDW+ 
SDW- 



91+92 < 0, 
-.gi +92 <0, 
-2gi + g2 - 93 > 0, 
-2.91 + .92 + 33 > 0, 
.92 + 53 > 0, 
92 - 53 > 0. 



(4.21a) 
(4.21b) 
(4.21c) 
(4.21d) 
(4.21e) 
(4.21f) 



The difference between the ladder and the parquet approximations shows up when there are more than one appre- 
ciable bricks in the problem. Then, the vertex /ssc contains not only the brick Cssci but other bricks as well, so Eqs. 
( 4.18 ) get coupled. This is the case, for example, for the ID spinlcss electrons, where the bricks C and Z are equally 
big, so they cancel each other in 7 (see Sec. 
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V. RESULTS OF NUMERICAL CALCULATIONS 



The numerical procedure consists of three consecutive steps; each of the m in volves solving differential equations by 
the fourth-order Runge-Kutta method. First, we solve parquet equations (4^) and (O) for the interaction vertices, 
which are closed equations. Then, we fi nd the triangu lar v ertices 7^, whose equations (4.7) and ( 4.12| ) involve the 
interaction vertices ji through Eqs. (4.8), (4.13), and (4.14). Finally, we calculate the susceptibilities Xi from Eqs. 
( 4.16 ), which depend on the triangular vertices 71. 

We select the initial conditions for the interaction vertices to be independent of the transverse momenta K.: 7i(/C, 
0) = gi. The momentum-independent interaction naturally appears in the Hubbard model, where the interaction is 
local in real space. In this C hap ter, the results are shown mostly for the repulsive Hubbard model without umklapp: 
9i = 92 = g, 53 = (Figs. ^^|l7|), or with umklapp: gi — 92 = gs = g (Figs. |l8|-[l^), where g is proportional the 
Hubbard interaction constant U. The absolute value of g (but not the sign of g) is not essential in our calculations, 
because it can be removed from the equations by redefining ^ to ^' = \g\£,. After the redefinition, we effectively have 
\g\ = 1 in the initial conditions. The actual value of \g\ matters only when the logarithmic variable ^' is converted 
into the temperature according to the formula T = ijLex'p{—2TTVF£,'/\g\). 

The initial independence of 7i(/C, C=0) on K, does not imply that this property is preserved upon renormalization. 
On the contrary, during renormalization, 7i(/C,^) develops a very strong dependence on IC and may even change sign 
in certain regions of the /C-space. We illustrate this statement in Fig. |l2| by showing typical dependences of 71 (/C, 
and 72(/C,^) on the average momentum py = (k^^ + k^^)/2 of the incoming electrons at ki = ks and k2 = k^ after 
some renormalization (^ = 1.4). In Figs. 12-141 the upper and lower limits on the horizontal axes are the boundaries 



±4°^ of the flat region on the Fermi surface, which are set to ±1 without loss of generality. One can observe in Fig. |1^ 
that the electron-electron interaction becomes negative (attractive) at large Py, even though initially it was repulsive 
everywhere. 



Mathematically, the dependence of 7i(/C,^) on /C arises because of the finite limits of integration, [—k'jj'\ 



imposed on the variables fcy'^-' and ky'"' in Eqs. (|4.3|) 



zero, fcy'^'' may change from —k. 



(0) 



while k: 



(B) 



.(1) 



to be confined to a narrower interval [—kf^ + 2py,k'^'] to ensure that ky"' = 2py — ky^' stays within [—k'^'^ky 
This difference in the integration range subsequently generates the dependence of 7i(AC, ^) on py and, more gen erally, 
on JC. Since many channels with different geometrical restrictions contribute to dji{K.,£,)/d^ in Eqs. ( [4.3| ), the 
resultant dependence of 7i(/C,^) on the four-dimensional vector IC is complicated and hard to visualize. Because of 
the strong dependence of 7i(/C,C) on /C, it is not possible to describe the 2D system by only three renormalizing 
charges 7i(C)i 72(0; ^^'^ 73 (Oj the ID case. Instead, it is absolutely necessary to consider an infinite number 

of the renormalizing charges "fi{IC,^) labeled by the continuous variable IC. This important difference was neglected 
in Ref. ||^, where the continuous variable IC was omitted. 

Having calculated 7i(A^,0, we solve Eqs. (4.7) and (4.12) for the triangular vertices (the order parameters) 
T{ky,qy,S,), which depend on both the relative {ky) and the total (qy) transverse momenta. We find numerically 
that the order parameters with qy = diverge faster than those with qy ^ 0. This is a natural consequence of the 
integration range restrictions discussed above. For this reason, we discuss below only the order parameters with zero 
total momentum qy = 0. Wc select the initial conditions for the symmetric and antisymmetric order parameters in 
the form: 



For example, in Eq. (4.3a), when py = (fcj 
stays in the same interval. However, when py 



,(0)l 



JB) 



(2) 

ky ')/2 equals 
> 0, has 

^(0) 1,(0)1 



7-+(fc„e=0) = l, 7^-(fc„e=0) 



ky. 



(5.1) 



In Figs. [T^ and |l^, we present typical dependences of the superconducting and density-wave order parameters on 
the relative momentum ky at the same renormalization "time" f = 1.4 as in Fi g, The singlet antisymmetric 
component (^g^) dominates among the superconducting order parameters (Fig. O), whereas the symmetric SDW 
order parameter (7^^^) is the highest in the density- wave channel (Fig. p^ . 

Having calculated the triangular vertices T, we find the susceptibilities from Eq. (4.16). The results are shown 
in Fig. |l^. The symmetric SDW has the fastest growing susceptibility Xsdw which diverges at ^sdw = 1-76. 
This divergence indicates that a phase transition from the metallic to the antiferromagnetic state takes place at 
the transition temperature TgDw = Mexp(— 27rz)p'^sDw/5)- A similar result was obtained in Ref. [ p9| by analyzing 
the convergence radius of the parquet series in powers of g ^. In the ladder approximation, the SDW instability 
would take place at CgDW ~ ^/92 — 1, as follows from Eqs. ( 4.14 ) and ( 4.18c ). Since ^sdw > ^sdW' ^^^'^ transition 



temperature TsdWi calculated in the parquet approximation, is lower than the temperature TsdW' calculated in the 
ladder approximation: Tsdw < ^smv- "^^^ parquet temperature is lower, because competing superconducting and 
density-wave instabilities partially suppress each other. 
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Thus far, we considered the model with ideally flat regions on the Fermi surface. Suppose now that these regions 
are only approximately flat. That is, they can be treated as being flat for the energies higher than a certain value 
^-cutoff, but a curvature or a corrugation of the Fermi surface becomes appreciable at the smaller energies E < i?cutoff- 
Because of the curvature, the Fermi surface does not have nesting for E < £^cutofi; thus, the density- wave bricks in 



the parquet equations (4.3) stop to renormalize. Formally, this effe ct can be taken into account by introducing a 



cutoff ^cutoff — (l/27rui?) ln(/x/i?cutoff), so that th e r.h .s. of Eqs. (4.3c )-(4.31) for the density-wave bricks are replaced 



by zeros at ^ > ^cutoff- At the same time, Eqs. (4.3a) and (4.3b) for the superconducting bricks remain unchanged, 
because the curvature of the Fermi surface does not affect the superconducting instability with qy = 0. The change of 
the renormalization equations at ^cutoff is not a completely rigorous way |Q to take into account the Fermi surface 
curvature; however, this procedure permits obtaining explicit results and has a certain qualitative appeal. For a more 
rigorous treatment of the corrugated Fermi surface problem see Ref. |3^] . 

In Fig. H^, we show the susceptibilities calculated using the cutoff procedure with ^cutoff = 1-4. The density- wave 
susceptibilities remain constant at ^ > ^cutoff- At the same time, Xssc(0 diverges at Cssc ~ ^-^^ indicating a 
transition into the singlet superconducting state of the d-wave type. Thus, if the SDW instability is suppressed, the 
system is unstable against formation of the d-wave superconductivity. This result is in agreement with the conclusions 
of Refs. ||,|^,||. 

From our numerical results, we deduce that the dependence of C^g^ on ^cutoff is linear: ^gg^ — a — &Ccutoff with 
b = 2.06, as shown in the inset to Fig. Converting ^ into energy in this relation, we flnd a power law dependence: 

Tssc'^T^r—- (5-2) 

^cutoff 



Eq. (5.2) demonstrates that increasing the cutoff energy i^cutoff decreases the temperature of the superconducting 
transition, Tgg(^. Such a relation can be qualitatively understood in the following way. There is no bare interaction in 
the superconducting d-wave channel in the Hubbard model, so the transition is impossible in the ladder approximation. 
The growth of the superconducting d-wave correlations is ind uced by the growth of the SDW correlations, because 



the two channels are coupled in the parquet equations ( 4.18 ). If ii^cutoff is high, the SDW correlations do not have 



enough renormalization-group "time" ^ to develop themselves because of the early cutoff of the density- wave channels; 
thus, Hence, decreasing -Ecutoff increases Tg^^-,. However, when E'cutoff becomes lower than TsdW: the 

SDW instability overtakes the superconducting one. Corresponding phase diagram is shown in Fig. |l^. Generally 
speaking, the phase diagram plotted in the energy variables, as opposed to the logarithmic variables ^, may depend 
on the absolute value of the bare interaction constant \g\. In Fig. |l^, we placed the points for the several values of 
g = 0.3, 0.4, and 0.5; the phase boundary does not depend much on the choice of g. The phase diagram of Fig. |l^ 
qualitatively resembles the experimental one for the high- Tc superconductors, where transitions between the metallic, 
antiferromagnetic, and superconducting states are observed. The value of -Ecutoff may be related to the doping level, 
which controls the shape of the Fermi surface. Taking into account the crudeness of our approximations, detailed 
agreement with the experiment should not be expected. 

We perform the same calculations also for the Hubbard model with umklapp scattering {gi — g2 — 93 — As 
one can see in Fig. ^ where the susceptibilities are shown, the umklapp process does not modify the qualitative 
picture. The leading instability remains the SDW of the symmetric type, which is now also symmetric with respect 
to the umklapp scattering, whereas the next leading instability is the singlet d-wave superconductivity. The SDW 
has a phase transition at ^jnwo- = 0.54, which is close to the ladder result fsDW+ ~ 1/(32 + gs) = 0.5. Some of 
the susceptibilities in Fig. Iia coincide exactly, which is a consequence of a special SU(2)xSU(2) symmetry of the 
Hubbard model at the half filling |^^. The phase diagram with the energy cutoff (Fig. |l^) is similar to the one 
without umklapp (Fig. p7| ), but the presence of the umklapp scattering decreases the transition temperature of the 
d-wave superconductivity. 

An important issue in the study of the ID electron gas is the so-called g-ology phase diagram, which was constructed 
for the flrst time by Dzyaloshinskii and Larkin |l^]. They found that, in some regions of the (51,(72,53) space, the 
ID electron system develops a charge or spin gap, which is indicated by divergence of 7i(^) with increasing ^. In the 
region where none of the gaps develops, the Luttinger liquid exists. It is interesting whether such a region may exist 
in our 2D model. To study the phase diagram of the 2D system, we repeat the calculations, systematically changing 
relative values of gi, 52, and 53. From the physical point of view, the relative difference of 51, 52, and 53 roughly 
mimics dependence of the interaction vertex on the momentum transfer. As an example, we show the susceptibilities 
in the case where gi = 2, 52 = 1, and 53 = in Fig. In this case, the leading instabilities are simultaneously the 
triplet superconductivity of the symmetric type (TSC+) and the spin-density wave. 

For all studied sets of gi, we find that the leading instabilities calculated in the parquet and the ladder approxi- 
mations always coincide. (We do not introduce the energy cutoff here.) Thus, the parquet effects do not modify the 
g-ology phase diagram of the 2D model derived in the ladder approximation, even though the transition temperatures 
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in the parquet approximation are always lower than those obtained in the ladder approximation. In that sense, the 
parquet corrections are much less important in the 2D case than in the ID case. From the mathematical point of 
view, this happens because a leading divergent brick develops a strong dependence on the transverse momenta /C and 
acquires the so-called mobile pole structure po| , ^ , ^ : 

^(/COoc^^. (5.3) 



The name "mobile pole" is given, because the position of the pole in ^ in Eq. (5.3), ^c(^), strongly depends on the 
momenta /C. It was shown in Rcfs. H^,^,^ that, because of the mobility of the pole, the leading channel decouples 
from the other channels, and the parquet description effectively reduces to the ladder one, as described at the end of 



Sec. [V. The phase diagram of the 2D system in the ladder approximation is given by Eqs. ( 4.21 ). It follows from 



Eqs. (4.21) that every point in the (31,52,53) space has some sort of instability. Thus, the Luttingcr liquid, defined 
as a nontrivial metallic ground state where different instabilities mutually cancel each other, does not exist in the 2D 
model. 

Generally speaking, other models may have different types of solutions of the fast parquet equations, such as 



immobile poles P0|] or a self-similar solution |23|, the latter indeed describing some sort of a Luttinger liquid. In our 
study of a 2D model with the van Hove singularities | |3^ , we found a region in the 5-space without instabilities, where 
the Luttinger liquid may exist However, we find only the mobile-pole solutions in the present 2D model. 

VI. CONCLUSIONS 

In this paper we derive and numerically solve the parquet equations for the 2D electron gas whose Fermi surface 
contains flat regions. The model is a natural generalization of the ID electron gas model, where the Luttinger liquid 
is known to exist. We find that, because of the finite size of the flat regions, the 2D parquet equations always develop 
the mobile pole solutions, where the leading instability effectively decouples from the other channels. Thus, a ladder 
approximation is qualitatively (but not necessarily quantitatively) correct for the 2D model, in contrast to the ID case. 
Whatever the values of the bare interaction constants are, the 2D system always develops some sort of instability. 
Thus, the Luttinger liquid, defined as a nontrivial metallic ground state where different instabilities mutually cancel 
each other, does not exist in the 2D model, contrary to the conclusions of Refs. ]2^ , |27t . 

In the case of the repulsive Hubbard model, the leading instability is the SDW, i.e., antiferromagnetism If 
the nesting of the Fermi surface is not perfect, the SDW correlations do not develop into a phase transition, and the 
singlet superconductivity of the d-wave type appears in the system instead. These results may be relevant for the 
high- Tc superconductors and are in qualitative agreement with the findings of Refs. |25 3C| , |3^ . 

In the bosonization procedure |4|-[lll, a higher-dimensional Fermi surface is treated as a collection of flat patches. 
Since the results of our paper do not depend qualitatively on the size of the flat regions on the Fermi surface, the 
results may be applicable, to some extent, to the patches as well. Precise relation is hard to establish because of the 
infinitesimal size of the patches, their different orientations, and uncertainties of connections between them. On the 
other hand, the bosonization procedure seems to be even better applicable to a flat Fermi surface, which consists of 
a few big patches. Mattis and Hlubina |27| followed that logic and claimed that the flat Fermi surface model is 
exactly solvable by the bosonization and represents a Luttinger liquid. The discrepancy between this claim and the 
results our paper indicates that some conditions must restrict the validity of the bosonization approximations. Luther 
gave a more sophisticated treatment to the flat Fermi surface problem by mapping it onto multiple quantum chains 
p8[ . He found that the bosonization converts the interaction between electrons into the two types of terms, roughly 
corresponding to the two terms of the sine-Gordon model: the "harmonic" terms {dip/dx)'^ and the "exponential" 
terms exp(i(/3), where is a bosonization phase. The harmonic terms can be readily diagonalized, but the exponential 
terms require a consistent renormalization-group treatment. If the renormalization-group equations were derived in 



the bosonization scheme of |28|, they would be the same as the parquet equations written in our paper, because 
the renormalization-group equations do not depend on whether the boson or fcrmion representation is used in their 
derivation [|||. 

Long time ago, Luther bosonized noninteracting electrons on a curved Fermi surface ||4^ ; however, the interaction 
between the electrons remained intractable because of the exponential terms. The recent bosonization in higher 
dimensions managed to reformulate the problem in the harmonic terms only. This is certainly sufficient to 

reproduces the Landau description of sound excitations in a Fermi liquid Q]; however, it may be not sufficient to 
derive the electron correlation functions. The validity of the harmonic approximation is hard to trace for a curved 
Fermi surface, but considerable experience has been accumulated for the flat Fermi surface models. 
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In the model of multiple ID chains without single-electron tunneling between the chains and with forward scattering 
between different chains, the bosonization produces the harmonic terms only, thus the model can be solved exactly 
1^, 45|. However, a slight modification of the model by introducing backward scattering between different chains 
1 2^, 16 1 or interaction between four different chains [|T| adds the exponential terms, which destroy the exact solvability 
and typically lead to a CDW or SDW instability. Even if no instability occurs, as in the model of electrons in a 
high magnetic field j2^, the fast parquet method shows that the electron correlation functions have a complicated, 
nonpower structure, which is impossible to obtain within the harmonic bosonization. Further comparison of the fast 
parquet method and the bosonization in higher dimensions might help to establish the conditions of applicability of 
the two complementary methods. 

The work at Maryland was partially supported by the NSF under Grant DMR-9417451, by the Alfred P. Sloan 
Foundation, and by the David and Lucile Packard Foundation. 
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FIG. 1. Dispersion law of ID electrons. The states in the shaded range of the momentum kx are occupied by electrons. 
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FIG. 2. Bare superconducting and density-wave susceptibilities. The solid and dashed lines represent the Green functions of 
the + and — electrons. The wavy lines represent incoming momentum and energy. 
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FIG. 3. Some parquet corrections to the vertex of interaction between electrons, 7, which is shown as a circle. The dots 
represent the bare interaction vertex g. The expressions beneath the diagrams represent the values of the corresponding 
diagrams. 
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FIG. 4. Decomposition of the interaction vertex 7, shown as a circle, into superconducting and density-wave bricks, shown 

as rectangles, in the spinless case. 
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FIG. 8. Vertices of interaction between electrons with spin: backward (71), forward (72), and umklapp (73) scattering. The 
thin solid lines inside the circles indicate how spin is conserved. 
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FIG. 10. Parquet equations for the bricks for electrons with spin. The variables ki and k2 (fca and fc4) represent momenta 
of incoming (outgoing) electrons, whereas the variables kA and ka represent intermediate momenta that should be integrated 
over. 
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FIG. 11. Fermi surface of a 2D electron gas. The thick hnes indicate flat regions on the Fermi surface. 
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FIG. 12. Interaction vertices 71 (fcy^' , k^y'^ ; A^'^' , k^f^ ; ^) and 72 (/cy^' , fc^^^ ; k'^y^ , k^y^ ; ^) as functions of the average momentum 
Pv = {ky^ + k^y^)/2 of the incoming electrons at k^y^ = ky^\ ky^^ = and { = 1.4. 
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FIG. 13. Superconducting order parameters Tsc{ky,qy,i) as functions of relative momentum ky at = and ^ = 1.4. 
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FIG. 14. Density-wave order parameters Tow{ky,qy,^) as the functions of relative momentum ky at qy = and ^ = 1.4. 
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FIG. 15. Evolution of susceptibilities Xi(^) in the repulsive Hubbard model without umklapp. 
i = CsDw = 1.76. 



XsDw(0 diverges at 
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FIG. 16. Evolution of the four leading susceptibilities in the Hubbard model without umklapp, with the cutoff of all den- 
sity-wave channels at ^ > ^cutoff = 1-4. The susceptibility of antisymmetric singlet superconductivity, Xssc diverges at 
^ = esse = 2-44. 
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FIG. 17. Phase diagram of the Hubbard model without umklapp illustrates the dependence of the d-wave superconducting 
transition temperature on the cutoff energy Scutoff for different g. The inset shows the same dependence in the logarithmic 
variables ^ggQ and ^cutoff. The solid line in the inset is a fit: SJ^qq = a — b ^cutoff, a = 5.33, b = 2.06. 
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FIG. 18. Evolution of generalized susceptibilities x<($) in the Hubbard model with umklapp scattering. Xsdw+(^) diverges 
at C = isDw+ = 0-54. 
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FIG. 19. Phase diagram of the Hubbard model with umklapp scattering illustrates the dependence of the critical temperature 

^ssc cutoff energy -C/cuto 

£F for different g. The inset shows the dependence of ^ggQ on ^cutoff- The solid line in the inset 



is a fit: ^gsc 



b Ccutoff, a = 4.5, b = 7.32. 
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FIG. 20. Evolution of susceptibilities in the case where gi = 2, g2 = 1, and = 0. Xsdw(^) Xtsc(^) diverge at ^ = 3.48. 
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